Functional renormalization group study 
of the interacting resonant level model in and out of equilibrium 
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We investigate equilibrium and steady-state non-equilibrium transport properties of a spinless 
resonant level locally coupled to two conduction bands of width ~ Y via a Coulomb interaction U 
and a hybridization t' . In order to study the effects of finite bias voltages beyond linear response, 
a generalization of the functional renormalization group to Keldysh frequency space is employed. 
Being mostly unexplored in the context of quantum impurity systems out of equilibrium, we bench- 
mark this method against recently-published time-dependent density matrix renormalization group 
data. We thoroughly investigate the scaling limit T — > oo characterized by the appearance of power 
laws. Most importantly, at the particle-hole symmetric point the steady-state current decays like 
J ~ V~ aj as a function of the bias voltage V S> t' , with an exponent aj(U) that we calculate to 
leading order in the Coulomb interaction strength. In contrast, we do not observe a pure power-law 
(but more complex) current-voltage-relation if the energy e of the resonant level is pinned close to 
either one of the chemical potentials ±V/2. 

PACS numbers: 71.10Pm, 73.63Kv 



I. INTRODUCTION 

Experiments on nanostructures represent a highly- 
active field of research. Whereas transport properties 
can be measured straightforward beyond linear response, 
a theoretical approach to quantum impurities out of equi- 
librium is challenging in presence of Coulomb interactions 
which are ubiquitous in low-dimensional systems. Over 
the last years, a great variety of both numerical as well as 
analytical methods was developed to study correlation ef- 
fects on the non-equilibrium dynamics of or steady-state 
current through quantum dots. Ranking among those 
are exact Bethe ansatz solutions,— £ perturbative renor- 
malization group schemes^— quantum Monte Carlo^ 
a real-time path integral approach^ Hamiltonian flow 
equations^ as well as the time-dependent numerical 
(NRG) 11 ' 12 and density matrix renormalization group 
(DMRG)i£~— frameworks. While all these methods had 
long-standing success in computing linear-response prop- 
erties of quantum impurity systems, the non-equilibrium 
situation is still a newly-emerging and thus rather unex- 
plored field. 

The functional renormalization group (FRG) provides 
an a priori exact re-formulation of a correlated many- 
particle problem in terms of coupled flow equations for 
irreducible vertex functions of arbitrary ordcr^i 7 . In the 
context of quantum dots in equilibrium, even a very 
simple way to truncate this infinite hierarchy, which 
can be viewed as a kind of RG enhanced Hartree-Fock 
approximation, allows for accurately describing the ef- 
fects of small to intermediate (and sometimes even large) 
Coulomb interactions U very flexibly and with minor nu- 
merical effort. Most importantly, the zero-temperature 



linear conductance has been computed in good agree- 
ment with NRG reference data for a variety of quantum 
dot geometries. 18 ' 19 Employing a more elaborate trunca- 
tion scheme (where one accounts for the frequency de- 
pendence of the two-particle vertex) allows for calculat- 
ing finite-energy properties like the density of states at 
least for intermediate values of U (for the single impurity 
Anderson model) ] 20 ' 21 

In contrast to the regime of linear transport where 
the strength and limitations of the functional RG in 
the context of quantum impurity systems have been ex- 
tensively investigated, there are only few works on the 
case of non-equilibrium. Even though it was proven 
possible in the steady-state limit to derive an infinite 
hierarchy of flow equations in Keldysh frequency space 
which are structurally identical to those on the Matsub- 
ara axis^ 2 - - — little is known on how different approxi- 
mation (i.e., truncation) schemes succeed or fail to de- 
scribe out-of-cquilibrium physics of correlated quantum 
dots. The first aim of this paper is to partly fill this gap 
by benchmarking functional RG calculations for a very 
simple impurity model (namely the interacting resonant 
level model) against recently-published linear-response 
and time-dependent DMRG data^^ 

Quite generally, the interacting resonant level model 
(IRLM) describes a single localized level (with an energy 
e) coupled to a bath of delocalized states (featuring a 
bandwidth ~ T) both by a local Coulomb repulsion U 
and a hopping matrix element t'. It was initially intro- 
duced four decades ago to study the equilibrium physics 
of mixed- valence compounds, and observables were com- 
puted by mapping to the anisotropic Kondo model (and 
using results available for the latter) or by pcrturba- 
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FIG. 1: (Color online) Schematic presentation of the two- 
channel interacting resonant level model studied in this paper. 



tive RG calculations.— ~— The two-channel version of the 
IRLM, which has gained considerable interest within the 
past few years , 2 ' 26 i 32 ~ — represents a very simple impu- 
rity model to describe charge fluctuations and investi- 
gate non-equilibrium transport (driven by a bias volt- 
age V = (hl + A*ii)/2 between the two baths) through 
a quantum dot. Most notably, accurate time-dependent 
DMRG data was recently published by Boulat et al. (for 
fairly large values of t'/T) and provides the aforemen- 
tioned benchmark for the functional RGj 2 ^ The oppo- 
site (so-called scaling) limit of the bandwidth T being 
much larger than all other energy scales was charac- 
terized by the appearance of universal power laws (by 
approximate approaches each having its advantages and 
shortcomings) i 26 ' 34 i 35 In particular, the current through 
the system was found to decay like J ~ V~ aj ^ u ' c ' for 
r ^> V ^s> t', both for the impurity energy being small 
(e ^ I/) 34 i 35 or close to one of the chemical potentials 
(e ~ Hlm, aj(U, e /j, L>R ) = atj(U, e = 0)/2)M Having 
explored its own strength and drawbacks in comparison 
with the DMRG reference as well as with new equilibrium 
NRG data, the functional renormalization group allows 
for systematically studying the scaling limit of the mi- 
croscopic IRLM for small to intermediate values of the 
Coulomb interaction U . It is the second aim of this paper 
to provide a consistent picture of the zero-temperature 
physics in this parameter regime from the FRG point of 
view, particularly in relation with prior results. 

This paper is organized as follows. We introduce the 
interacting resonant level model as well as the functional 
renormalization group approach in non-equilibrium in 
Sees. |ll] and IIII1 respectively. Sec. IIVI is devoted to the 
comparison of FRG results with time-dependent DMRG 
data. We systematically investigate the scaling limit in 
Sec.fVland conclude with a brief summary. A more elab- 
orate FRG approximation scheme is briefly discussed in 
the Appendix. 



II. THE MODEL 

The interacting resonant level model is depicted 
schematically in Fig. [TJ It describes a single spinless level 
of energy e as well as two (left and right) bath of delo- 
calized states: 

H imp = (e-U/2)4d 2 , H^ R = Y.vl c sk- (1) 



We model a local Coulomb interaction U and hopping t' 
between both parts by adding two distinguished neigh- 



boring sites, 
H v = U (d t 2 d 2 d\d 1 + d\d 2 d\d^j 

-t' (4^1 + 4 d 3 + H-C-) - U/2 (d^i + > 



(2) 



iJcoup 7= { d l C Lk + 4 c Rk + H - C - ) ' 



which are coupled to the bath of size N via 
t 



The characteristic energy scale (i.e., the bandwidth) of 
the latter is determined by the hybridization 

T = 7rp bath (w = 0)t 2 , (4) 

with pbath(w) being the local density of states. In order 
to explicitly compare with DMRG results, we model the 
baths as semi-infinite tight-binding chains with a nearest- 
neighbor hopping amplitude t and correspondingly 

PbathH = tt^V 7 ^ 2 - lo 2 9(2* - M) , T = t. (5) 

Z7TI 

The associated retarted Green function at the end of the 
isolated chain is given by 



ret 1 \uj- sgnHVw 2 - 4i 2 M > 2t 



As an alternative, one frequently employs completely 
structureless (wide-band) leads featuring a constant local 
density of states and 

flbathM = -iTPbath ■ (7) 

In our case, such realization is used to investigate the 
scaling limit T — > oo where details of the dispersion 
do not play any role (which one can can show explicitly 
within the FRG framework; see Sec. fV|) . 

The equilibrium statistics of the interacting resonant 
level model is determined by the usual grand canonical 
density operator p = exp(—[3H) featuring an inverse tem- 
perature /3 and equal chemical potentials /.jl = Pr = 0. 
The non-equilibrium situation is modeled by an initially 
separated system (H coup = 0) in a thermal bulk state 



(8) 



where pi mp denotes the density matrix of the isolated in- 
teracting three-site region, which we choose to be that of 
a vacuum configuration. At some time to, the coupling is 
switched on, and the time evolution for t > to is governed 
by the full Hamiltonian 



H = Ht 



imp 



Hu + H 



coup 



H 



bath 



H 



bath 



(9) 

In presence of a finite bias voltage V = 2fiL = — 2/i«, 
one does in general expect the system to relax to a non- 
thermal steady state independent of /5i mp at t — > oo, and 
this scenario is supported by time-dependent DMRG cal- 
culations for the problem at hand. 26 In this paper, we 
will focus exclusively on studying the steady state of the 
IRLM in the zero-temperature limit. 
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III. THE METHOD 

A. General idea of the functional RG 

The functional renormalization group implements Wil- 
son's general RG idea in terms of an infinite hierarchy of 
differential flow equations for single-particle irreducible 
vertex functions (such as the self-energy), which alto- 
gether represents an exact reformulation of the underly- 
ing many-particle problemJi The hierarchy is set up by 
introducing an infrared cutoff A into the non-interacting 
Green function Go, 

G (l , ;l)->Gft(l , ;l), Gq = °° = , G^ = G , (10) 

where the arguments are a shorthand for single-particle 
quantum numbers as well as either Matsubara frequen- 
cies (in equilibrium) or real frequencies and Keldysh in- 
dices (in non-equilibrium). Under the assumption of the 
existence of a steady state for the latter case (which al- 
lows for associating energies with time differences), such 
an infrared (energy) cutoff can be devised straightfoward, 
but this will be postponed to the next Section since its 
actual form is irrelevant for the time being. By virtue of 
the replacement (|10[) , every vertex function acquires a A- 
dependence, and both in and out of equilibrium one can 
derive the structurally same set of functional RG flow 
equations by straight-forwardly differentiating with re- 
spect to the cutoff parameter A. This can be technically 
achieved, e.g., by considering generating functionals, and 
the flow of the self-energy is given b y 18 i 22 ' 23 

3 A S A (l'; 1) = - £ 5 A (2; 2') 72 A (1', 2'; 1, 2) , (11) 

22' 

with the so-called single-scale propagator 

S A = -G A d A [Go] _1 G A = -d A G A , (12) 

and G A being the full Green function at scale A. Both the 
two-particle vertex 7^ and £ A itself enter the right-hand 
side of the differential equation (fTTj). Similarly, the flow 
of the 71-th order function 7 A is in general determined 
by all vertices J A <n - Integrating this infinite set of cou- 
pled differential equations from A = 00 (where all energy 
scales are suppressed and the many-particle problem be- 
comes trivial) down to A = (where one recovers the 
full energy spectrum) yields an in principle exact expres- 
sion for the self-energy of the system under consideration. 
In practice, however, one needs to devise a truncation 
scheme. In the main part of this paper, we focus solely 
on the flow equation pip with the two-particle vertex set 
to its initial value (j A=oc ~ U) , rendering £ a frequency- 
independent (effective) quantity. This Hartree-Fock-like 
approximation is correct at least to leading order in U 
and thus a priori justified in the limit of small Coulomb 
interactions. In equilibrium, it was proven to give reli- 
able results up to intermediate values of U (and even to 



capture aspects of Kondo physics) for the linear conduc- 
tance of a variety of quantum dot geometries! 18 ' 19 A more 
elaborate truncation schem e 20 ' 21 which accounts for the 
flow of the frequency-dependent two-particle vertex will 
be discussed briefly in the Appendix. 



B. Green functions &: Dyson equation 

1. Green functions in and out of equilibrium 

Linear-response transport properties of the interacting 
resonant level model can be computed conveniently using 
Matsubara Green functions 

G£M (r r a 2 (r)ot(0)) dr 

= [G^M- 1 -s^M]-. 1 , 

where o^(t) is a fcrmionic annihilation operator in the 
Heisenberg picture, T T denotes ordering with respect to 
the imaginary time t, and the self-energy E oq is associ- 
ated with the Coulomb interaction. In order to employ 
standard diagrammatic techniques in non-equilibrium 
(such as the very notion of vertex functions), we define 
long-time Green functions as 

Gy(w) = -i I e iuJt lim (t^ (iW(0)\ dt 

J to-t-oo \ ■> I p 

(14) 

with T c being the order operator on the Keldysh contour 
whose branches are characterized by indices a = ±. It is 
often more convenient to directly exploit causality (par- 
ticularly in approximate schemes which naturally con- 
serve this symmetry; see Sec. IIII C[) and work in the re- 
tarded, advanced, and Kelydsh basis: 

G rct = q— _ q-+ = ( G adv)t ^ G K = q-+ + q+- ^ 
E rct = s — + 5,-+ = ( E adv)t ^ = j,— + s ++ _ 

(15) 

In the next Section, we will explicitly derive the (non- 
interacting) Green functions of the IRLM in and out of 
equilibrium. 



2. Dyson equation 

For the problem at hand, the flow of the self-energy 
is determined by the Green functions of the interact- 
ing three-site system only. Using standard projection 
techniques^ the latter can be expressed in terms of a 
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finite matrix Dyson equation 
G rot H _1 = G^(w)- 1 - E rct (w) 



cquil. 



G eq (iw^ u + iOy 1 



G K H = G ret H[£ K H + £L h M]G ad >), 

(16) 

where the retarded Green function of the isolated impu- 
rity region is given by 

fuj + iO -t' \ 

flSfpM _1 = u-e + M -f , (17) 

\ -t' ui + iOJ 

and the self-energy associated with the bath reads 



1 





S bathM = -2i7rp bath t 2 I 



(18) 



1 - 



The latter are initially in a thermal state (for which the 
dissipation-fluctuation theorem holds) and described by 
the Fermi functions 



e {/3-^L.n)u _|_ I 



(19) 



In the next Section, we will derive the functional renor- 
malization group flow equations in order to compute an 
approximation for the self-energies S ret and E K (out of 
equilibrium) or S cq (for linear-response) which in each 
case incorporate the effects of the Coulomb repulsion. Af- 
terwards, one can calculate the current flowing through 
the interacting region as well as the zero-tempcraturc 
equilibrium conductance using the formulas (in units of 
e 2 /h = l)& 



Js = 2^ 2 J Pbath [/. (G+- - + G~+] du 



(20) 



G = 4r 2 |G^ = 0)| 2 , 



with the Green function's single-particle index being i = 
1 or i — 3 for the current at the left and right interface 
(s = L,R), respectively. Within all FRG approximation 
schemes employed to study the problem at hand, current 
conservation Jl = — Jr holds (whereas other symmetries 
may be violated; see Sec. MI Cl) . 



C. Flow equations 

In this paper, we will almost exclusively focus on con- 
sidering the flow of the self-energy only (i.e., truncat- 
ing the infinite RG hierarchy to leading order). This is 



achieved by setting the two-particle vertex 7 A to its ini- 
tial value in Eq. ([lip , which in equilibrium is nothing but 
the bare frequency-independent Coulomb interaction 

7 A (l', 2'; 1, 2) = ±U(3- 1 S(iuj v + ioj 2 , - iuJi - iuj 2 ) (21) 

for all permutations of the nearest-neighbor single- 
particle indices. An additional factor of imaginary i ap- 
pears in a real frequency representation, and the two- 
particle vertex is only non- vanishing if all Kclydsh indices 
a are equals 



a.2' 



Oi\ = a 2 



7 2 A (1', 2'; 1, 2) = ±aiUS(u> v + uo 2 , - cj x - u 2 ) . (22) 

As in Eq. (|21[) , different signs refer to symmetric or an- 
tisymmetric ways of ordering nearest-neighbor quantum 
numbers. As mentioned above, the self-energy obtained 
from this truncation scheme is frequency-independent 
(i.e., Hartree-Fock-like) and at least correct to first or- 
der in U. It contains, however, certain classes of higher- 
order contributions due to the underlying RG proce- 
dure and was shown to be reliable for small to inter- 
mediate Coulomb interactions and even to capture non- 
perturbative aspects of strong electronic correlations in 
the context of equilibrium quantum dot systemsi 18 ' 19 A 
more elaborate approximation scheme which accounts for 
the frequency-dependence of the two-particle verte x 20 i 21 
will be discussed briefly in the Appendix (for the IRLM 
in linear response). 



1. Equilibrium 

The last step in explicitly setting up the functional 
RG flow equations is to specify the form of the infrared 
cutoff. In equilibrium, low-energy degrees of freedom are 
most commonly suppressed by a sharp multiplicative O- 
function in Matsubara frequency space: 

G° q M -> Gg q ' A M = 9(M - A)Go q M , (23) 

and we will use such implementation throughout this pa- 
per in order to compute linear-response properties of the 
IRLM. The single-scale propagator S cq ' A , which solely 
determines the right-hand side of Eq. (TTT]) for a two- 
particle vertex set to its initial value, contains an at first 
sight ambiguous product of 0- and ^-functions. Evalu- 
ating the latter by means of Morris lemmaS yields 



5 oq ' A = 5(|w|-A) (Go q ) -£ cq 



= <5(M-A)G cq ' A , 
.(24) 

and the zero-temperature flow equations for the differ- 
ent independent self-energy components E^ q ' + t' = t' A , 



J 22 



e = e\, and £ 



cq,A 

11 



(introducing a notation 



where the interpretation as effective system parameters 



5 



becomes evident) read 
U 



fi»A*A = 
9a£A = 

9a<4 = 



Re 

7T 

-^Rc 

7T 

-^Re 

7T 



G^ A (zA) 



G^ A (*A) 



4^ = t' (25) 
= e (26) 
= (27) 
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The single-particle energy shift —U/2 appearing in the 
Hamiltonian of Eqs. ([TJ and © is cancelled by a contri- 
bution arising from analytically integrating from A = oo 
down to some arbitrarily large A — > oo. The ordi- 
nary coupled differential equations (f2"5]) - (j2"7) can be 
solved numerically (and sometimes even exactly; see 
Sec. |V} with minor effort, and we will discuss the linear- 
response physics described by this approximation scheme 
in Sees. El and M 



2. N on- equilibrium: sharp cutoff scheme 

One straightforward way of introducing a cutoff in non- 
equilibrium is the replacement analogous to Eq. 



d (u) -> G A (u) = 6(H - A) 



G+-(w) G++H 



(28) 

As before, the single-scale propagator S A which enters 
on the right-hand side of Eq. (fTTj) can be calculated by 
virtue of Morris lemmas 



*(M-A) 



G f 



S(\lu\ — A)G A , (29) 



and the flow of the self-energy is given by 



3 A S 



.A 

22 — 


LJ = ±A 


^ A H + G33-' A H 


.A 

11 — 


lU V G" 

27T ^ G 22 

LJ = ±A 


■ A (w) = 3 A S^' A 


.A 

H ~ 


iC/ V G" 

LJ = ±A 


■ A M, (ij) = (12,21,23,32) 


++.A _ 


-(a A s-- A ) t 




+ -.A = 


a A E- + - A = o , 





(30) 



complemented by the initial condition j^"* 00 = 0. 
Since non-equilibrium symmetry properties (particularly 
causality; see below) are not necessarily conserved by this 
approximation, it is not reasonable (but misleading) to 
interpret the different frequency-independent self-energy 
components as effective non-interacting system parame- 
ters. 

The very same sharp cutoff scheme defined by Eq. (f2"5)l 
has been previously applied to the single impurity An- 
derson model (SIAM)<2£ While qualitatively reproducing 



non-equilibrium features known, e.g., from perturbation 
theory, its major drawback (besides numerical problems) 
turned out to be the above-mentioned symmetry viola- 
tions. Most importantly, the causality condition 







(31) 



only holds to the truncation order (i.e., to leading U in 
the present context). Since the functional RG is a gen- 
erally non-conserving approximation^ the consequences 
of these violations of symmetries specifically associated 
with non-equilibrium are a priori unclear. Due to the 
lack both of reliable reference data as well as of an alter- 
native idea to introduce a conserving cutoff procedure, 
it was not possible to systematically address this ques- 
tion for the Anderson model at that point of time. Only 
recently^ Jakobs et al. introduced a scheme to suppress 
low-energy degrees of freedom (presented in the next Sec- 
tion) which does not violate causality in non-equilibrium 
(but features other shortcomings^), and thorough inves- 
tigations of the SIAM are under wayJi In this paper, 
we study the interacting resonant level model using both 
FRG cutoff schemes, particularly in comparison with ac- 
curate DMRG results^2£ 



3. N on- equilibrium: reservoir cutoff scheme 

The aforementioned alternative way to cut off low en- 
ergy modes within the functional renormalization group 
can be introduced on the Hamiltonian level as additional 
structureless reservoirs of zero chemical potential locally 
coupled to each site of the interacting region 



Hcnt — 



4EE44+H.C. 

V i=l k 



(32) 



where the hybridization A = wp C utt\ is used as the flow 
parameter. As we will show later on, this cutoff scheme 
preserves causality even after truncation, rendering it 
reasonable to directly work with retarded, advanced and 
Keldysh Green functions. The latter acquire a new self- 
energy-like term due to the additional reservoirs 

E£? t = -*A1 3 , S c K ut = -2iAtanh(/^/2)l 3 , (33) 

and the corresponding single-scale propagators of 
Eq. OH) read 

grot, A ^^rct,A^*rct,A ^gadv,A^t 

S K ' A = -d A [G 



rot. A /viK.A i V K i v K A ,^iadv,Al 



iret.A Ix^K.h , , v K \ f 

' ^bath ' ^cut U 



= s Ict - A (s 

Cr (2j -t- Zj bath + 2^ cut ) O 

^rrct, A/^»adv,A 



K "\ ^adv.A 
cut 

K 



+ 2itanh(/3o;/2)G rct ' A G a 



(34) 
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v/r v/r 

FIG. 2: (Color online) The steady-state current J as a function of the bias voltage V of the two-channel IRLM for large 
hoppings t' = 0.5r (in units of the bandwidth ~ T), zero impurity energy e, and various Coulomb interactions U. (a) 
Functional renormalization group results obtained from numerical integration of the self-energy flow equations (|3Up of the 
sharp cutoff scheme of Sec. IIII C 21 (b ) The same calculated with the reservoir cutoff approach of Sec. IIII C 31 Density-matrix 
renormalization group data of Ref. [2y for the same set of parameters is shown by symbols in the main part (for U = 0.3F only) 
as well as within the inset (where lines are a guide to the eye only). 



The zero-temperature flow of the effective system param- 
eters 



l 12 1 



i' = E 



/A 

23 

£a - e = E 

4 



^rot,A 
J 12 

rot. A 
23 

rot, A 
22 

iret,A 
J ll 



(ir-\adv,A 
^21 

v-»adv,A 
^32 

y*adv,A 
^22 

v»adv,A 
^11 



/ v ret,7 
^21 



^-»ret,A 
-^32 



*r->ret,A 

E 



22 

ret, A 
J ll 



yiret,A 
^33 



(35) 



can be derived straightforward by plugging the constant 
Coulomb interaction vertex given by Eq. (f2"2l into the 
general flow equation and rotating to the Keldysh 
basis. One obtains 



t 



/A- 



9a*23= I Sf^ A {uj) 

47T 



t 



/A— >oo 



23 



K,A/ 



9a4 = / S%i A (w) 



dui e\- 

4- 



t' (36) 

t' (37) 

e (38) 

0. (39) 



We note that E rct - A = £ adv ' A only holds within the first 
order FRG approach used in this paper. Another char- 
acteristic of the latter is that the Kelydsh component of 
the self-energy does not flow: 



K,A 



±g J [S^ A (w) + S adv < A H] duj 

\ iir — in = diag. components 
1 off-diag. components, 



(40) 



and the very same holds if one formally considers the 
'anti-causal' self-energy: 



d\ (E 



+ + ,A 



i .A 



KA 



— ,A\ 







(41) 



Thus, causality is not violated by this FRG approxima- 
tion scheme, providing the a posteriori justification to 
work in the basis of retarted, advanced, and Keldysh 
Green functions. The conservation of other symmetry 
properties (particularly of the current) follows straight- 
forwardly from the interpretation of the self-energy as 
effective (non-interacting) system parameters^ 



4-. From non- equilibrium to linear-response 

As mentioned above, we compute linear-response prop- 
erties of the IRLM using the Matsubara functional RG. 
The non-equilibrium formalism, however, is applicable 
for arbitrary bias voltages V and can thus be used to 
(approximately) describe equilibrium physics in the limit 

V — > 0. The latter is particularly simple within the reser- 
voir cutoff scheme introduced in Sec. IIII C 31 Namely, at 

V = the self-energy flow equation can be rewritten 
as (schematically omitting single-particle quantum num- 
bers) 

d A E ret ' A = t— [ d A G K - A (u)dco 
4?r J 

= ±^J tanh (^j d u [G^ A + G adv ' A ] du 
= ±— [G rct < A (0) + G adv ' A (0)l 



±-ReG cq < A (iA) 

7T 



(42) 



7 




'o 0.05 0.1 



e/r 

FIG. 3: (Color online) Linear-response conductance (in units 
of Go = e 2 /h) of the IRLM at t' = 0.1F as a function of 
the gate voltage e. The Figure shows a comparison between 
functional RG data obtained from the formalism of Sec. MI CTI 
(lines) and the DMRG results of Ref.[H (symbols). 



where we have used the fluctuation-dissipation theorem 
at zero temperature and assumed structureless wide- 
band leads with a frequency-independent self-energy 
given by ±iA, the latter allowing for the replacement 
d\ — > ±«<9 W . The resulting approximate flow equation 
(H2), however, is nothing but the Matsubara result at 
the same order of truncation.— In contrast, it is in gen- 
eral impossible to analytically analyze the flow within 
the sharp cutoff scheme in the linear-response limit sim- 
ilarly to Eq. (|42j) . and the same holds for the situation 
of baths featuring an energy-dependent local density of 
states. For the problem at hand, we have numerically 
studied the case of V = for those cases in order to en- 
sure that the results are always in agreement with the 
Matsubara formalism. 



IV. COMPARISON WITH DMRG 

In this Section, we show results both for the linear- 
response conductance and the steady-state current of the 
interacting resonant level model. We focus on the param- 
eter regime of large values of t' /T (particularly in non- 
equilibrium) for which reliable linear-response and time- 
dependent DMRG data was published recentlyj 25 ' 26 In 
order to explicitly compare with these results, we model 
our bath as infinite tight binding chains of bandwidth AT. 
The so-called scaling limit of T — > oo will be discussed 
extensively in Sec. [Vj 



A. Steady-state current 

The current J flowing between the two bath of the 
IRLM in presence of a finite bias voltage V is shown in 
Figs. (Hk) an d (b) for different values of the Coulomb in- 
teraction U, at the particle-hole symmetric point e = 0, 



and for fixed large hoppings t' = 0.5T. It was obtained 
from numerically integrating the non-equilibrium flow 
equations ([50)1 for the 9-approach (Fig. [U[a)) and (|56) 
- (|39| for the reservoir cutoff scheme (Fig. IUb)) as well 
the formula (f20|) . respectively. At U = 0, J increases lin- 
early for small bias voltages and saturates beyond some 
scale which is determined by t' (and will be quantified in 
Scc.fV]). In presence of a finite repulsive Coulomb interac- 
tion, it additionally features a regime V > t' 2 /V of nega- 
tive differential conductance (i.e., a current decreasing as 
the voltage is increased) which was frequently described 
in previous works (but still lacks a consistent physical 
explanation) i 2 ' 26 ' 34 ' 35 Since the hopping t' is fairly large, 
this decay is not governed by any specific (e.g., power- 
law-like) form. Most importantly, both FRG schemes 
show a satisfying agreement with DMRG data both for 
repulsive and attractive Coulomb interactions (see the 
symbols in the main part of Figs. EJa) and (b) as well 
as the inset to the latter), and the violation of causality 
prone to the sharp cutoff does not lead to unphysical re- 
sults. This indicates that even the most simple functional 
RG truncation scheme captures aspects of the essential 
non-equilibrium physics of the IRLM, giving confidence 
to use this approach in order to investigate the currcnt- 
voltage-relation more thoroughly in the so-called scal- 
ing limit r — > oo which cannot be accessed straightfor- 
ward within the time-dependent DMRG framework (see 
Sec. EJ). 



B. Linear response 

As a next step, we study the conductance G of the 
IRLM in the limit of linear response (see Fig. [3]) us- 
ing the equilibrium FRG introduced in Sec. IIII C II In 
the non-interacting case, its gate voltage dependence is 
given by a Lorentzian curve whose width is governed by 
t' (and shows power-law behavior if the latter becomes 
small compared to T; see Sec. |V|. In presence of small 
Coulomb interactions, the conductance is enhanced (and 
the Lorentzian widens) but eventually shrinks as U be- 
comes large. This effect was first observed in equilib- 
rium DMRG calculations^ and our approximate FRG 
data agrees quantitatively with those numerically exact 
results up to up (7/r w 1 and at least qualitatively for 
even larger interactions. 

In order to quantitatively compare with the DMRG 
data, we employed the IRLM with tight-binding leads 
featuring an energy-dependent local density of states. In 
that case, it is a priori impossible to analytically analyze 
the non-equilibrium FRG flow equations of Sec. IIII C 31 
in the linear-response limit, and the same holds for 
the sharp cutoff scheme in any case. We have there- 
fore checked numerically that the linear-response conduc- 
tance obtained from both non-equilibrium frameworks 
agrees with that of the Matsubara formalism (the latter 
being shown in Fig. [3]). 

All in all, the comparison with DMRG data establishes 
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u/r 

FIG. 4: (Color online) The exponents a.j and a x govern- 
ing the scaling-limit power-law behavior of the current (for 
t' <C V) and of the susceptibility (for V -C t'), respectively. 
To leading order, both quantities are given by a = 2U/ttV. 
For the equilibrium exponent a x , symbols show numerical 
renormalization group reference data. 



the very simple (Hartree-Fock-likc) FRG approach as a 
satisfying tool to investigate transport properties of the 
interacting resonant level model at small to intermediate 
Coulomb interactions in and out of equilibrium. 



V. SCALING LIMIT 

In this Section, we investigate the situation where the 
characteristic energy T of the bath of delocalizcd states 
(i.e., the bandwidth) is much larger than all other en- 
ergy scales. This so-called scaling limit was addressed 
in several prior works (with a special focus to non- 
equilibrium) and is supposed to be governed by universal 
power lawsi 34 i 35 We will particularly discuss the predic- 
tions of our FRG approximation scheme in relation with 
those earlier results. 

In the limit of large T, universal features (such as 
the exponents governing power laws) should only depend 
on the density of states at the chemical potential but 
not on the precise energetic structure of the bath. We 
have numerically checked that this indeed holds at least 
within our FRG approximation and will thus for simplic- 
ity model the leads as completely structureless through- 
out this Section. 



A. Renormalization of t' 

It is instructive to study the renormalization of the 
system parameters before discussing physical quantities. 
In the most simple case of equilibrium and zero impurity 
energy e = 0, the only remaining flow equation (|25p for 
of the hopping amplitude t' A takes the form 



/A 



u 



t ,A /r 



nr (A/r) 2 + A/r + 2(t' A /r) 2 



(43) 



with the initial condition given by t' A ^°° = t 1 . Thus, the 
flow of t lA starts below some ultraviolet scale A ~ r and 
is cut by itself in the infrared regime. More precisely, 
the differential equation (|4"3"]l can be solved analytically 
in the limit of small U/T and i'/F by setting t' A in the 
denominator to its initial value. The result reads 



t/A=0 



V 



(44) 



and the renormalizcd hopping scales as a power law with 
its bare value. Higher-order corrections to the exponent 
(which in contrast to the first order do not necessarily 
need to have the right prefactors as our truncated FRG 
is only correct to leading U) can be determined by nu- 
merically integrating the full Eq. (j43]l and fitting a power 
law to the solution t' A=0 (t'). As a side remark, we note 
that calculating the self-energy by first order perturba- 
tion theory in U (which can be most simply done by 
completely discarding the feedback of t' A on the right- 
hand side of Eq. (|4"3")l ) yields a logarithmic divergence in 
the bare system parameter^ 3 - 



j.1 pert 
t> 



ttF vr 



(45) 



indicating the necessity to employ an RG-based frame- 
work. 

Using similar analytical arguments at finite impurity 
energies and out of equilibrium, one can show that to 



leading order in U the flow equations for t a 3 
take the general formic 



Ki2.23l 2 /r 



A — a IV A + r A + r A 

dATl > 3 ~ ^^r(F/2 T e) 2 + (A 



(46) 



for A smaller than the ultraviolet cutoff ~ T. Thus, the 
effective hopping amplitudes ti j3 = t a ^° scale as a power 
law with each of the bare system parameters V, e, and t , 
provided that particular one is much larger (the meaning 
of which will be quantified in the next Section) than the 
other two but still much smaller than the bandwidth T: 



n,3 
7^ 



[(f'2)-S+°(f 2 ) V ,e < f <r 

v -S+o(c/ 2 ) t > }€ « y« r 

e -^+ ^ 2 ) V,t' « e «T 



(47) 



where higher-order corrections can again be extracted nu- 
merically. Such power-law behavior with respect to the 
voltage V and hopping t' was previously described using 
field-theoretical model a 26 ' 34 or perturbative renormaliza- 
tion group treatments^ and the corresponding expo- 
nents agree with the FRG result to leading order in U 
(despite some difficulties in relating the parameters of a 
continuum model to our microscopic ones). Moreover, 
it is suggested in Ref. HH that one can understand non- 
equilibrium properties from equilibrium only, which in 
the extreme limits mentioned above is supported by our 
observation that all bare system parameters (including 
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FIG. 5: (Color online) The current as a function of the voltage 
in the scaling limit F 4 oo obtained numerically from the 
sharp (solid lines) and reservoir (dashed lines) FRG cutoff 
schemes, respectively. Beyond some cross-over scale Tk, J 
decays as a power law of the voltage V over several orders 
of magnitude. The latter manifests as a constant logarithmic 
derivative dln(J/r)/dln(V/T). 



the voltage) can be interpreted in terms of a simple in- 
frared cutoff which automatically appear within the FRG 
framework (and do not have to be introduced by hand). 
In Sec. IV C[ we show that this altogether simple picture 
breaks down if the impurity energy is pinned to cither 
one of the chemical potentials (e s ±V/2). 



B. Susceptibility 

In this Section, we investigate how the power laws de- 
scribed above manifest in equilibrium observablcs. Two 
quantities associated with linear-response transport are 
the resonance width of the conductance G(e) and the 
charge susceptibility 



X 



d(n) 



dc 



-=o 



ttT k 



(48) 



with (n) being the average occupation of the resonant 
level. The latter can conveniently be computed within 
the equilibrium FRG formalism by integrating over the 
Matsubara Green function. The inverse susceptibility 
can be used to define a unique scale Tk which governs the 
low-energy linear-response physics4£ For small U, one 
can show that both the width of the conductance and \ 
are solely determined by the renormalizcd hopping t' A=0 
computed in the previous Section and thus governed by 
a power law in the linear-response limit V <C Tk' 



X 



-2a„ 



91 J 

,a x = — + 0(U 2 ) 7 (49) 



where the higher-order corrections to the exponent are in- 
fluenced by the flow of all self-energy components and can 
only be extracted numerically (by integration of Eqs. (|25[) 
- (1271) and subsequent power-law fitting; for the result 



see Fig. [4}. As mentioned above, one can in general not 
expect to obtain the right prefactor even of the second- 
order term since our truncated FRG scheme is correct 
only to leading U. However, going beyond first order 
allows both to specify the regime where the exponent is 
purely linear as well as for a quantitative comparison with 
other resultsi 34 ' 35 In our case, the equilibrium numerical 
rcnormalization group can be straightforward employed 
to compute a x with high accuracy, thus providing an 
additional benchmark for the functional RG in the scal- 
ing limit 4£ One observes that the very simple (Hartree- 
Fock-like) FRG approximation scheme of Eqs. (|2"5j) - (|2"T)) 
shows satisfying agreement with the NRG reference even 
for intermediate U/Y (sec Fig. 2]) 42 



C. Current 

The most interesting transport property of the IRLM 
in non-equilibrium is the current. One can show that for 
small U and large Y ^> t',V,e it is determined by the 
renormalizcd hoppings n,3, and an approximation to the 
latter can be derived from Eq. (|46[) . This yields 



J 



4:TxT 3 

n + r 3 



arctan 



fV/2-e 

\Tl + T 3 



arctan 



V/2 



Ti + T 3 



(50) 

In the following, we discuss this result as well as the nu- 
merically obtained current specifically for impurity ener- 
gies in the middle between or close to either one of the 
chemical potentials of the bath (e <C V and e « ±V/2, 
respectively). 



Zero impurity energy 



For large voltages V ^> Tk,£, Eq. 
power law 



T K \V J irY V ' 



50]) describes a 



(51) 



with an exponent a rthat to leading order agrees with 
the results of Refs. [H, [H, and HE Beyond the limit 
of small U, aj as well as the current itself can only be 
computed numerically by integrating the flow equations 
PP ]! and ([55 ]) - ([55 ]) . respectively. Within both FRG non- 
equilibrium cutoff schemes, we observe that in agreement 
with our analytics J generally features a linear increase 
crossing to a power-law decay (and thus a constant loga- 
rithmic derivative) at scale V ~ Tk^ but eventually sat- 
urates as one approaches the bandwidth Y (see Fig. [5]) . 
The numerically-determined exponent a,j is depicted in 
Fig. |4j It is purely linear up to sizable Coulomb interac- 
tions U £3 r, and our simple FRG approximation scheme 
does thus not contain higher-order corrections to the ex- 
ponent in case of non-equilibrium. 
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u/r 



0.2 



2.2 



2.22 



V/e 



FIG. 6: (Color online) The same as Fig. [5j but for a finite 
impurity energy e. The data was obtained by numerically 
integrating the flow equations (|36[) - <|39[) . The current is 
suppressed for small voltages V -C e but crosses over to a 
power-law decay when e is moved below the chemical poten- 
tial. Inset: The quantity a rcs governing the behavior close 
to the aforementioned resonance condition e ~ ±V/2 (see 
Eq.[52j. 



As a passing comment, we note that in order to actu- 
ally observe pure power-law behavior of the current, the 
voltage needs to be in a regime with Tk and 
the bare hopping amplitude t' typically has be chosen of 
the order of t' 2 /V 10" 6 (since T K ~ t' 2 /V at small U). 
By successively increasing t' we observe that the regime 
of voltages characterized by a power law shrinks until 
eventually for t' = 0.1T the logarithmic derivative only 
features a local minimum (with a value giving a rough es- 
timate of the exponent, though) close to F«r. Within 
the DMRG framework of Ref. [26|, such large hoppings 
were used for a power-law fit of the current. Even though 
it is certainly numerically demanding, treating smaller t' 
in non-equilibrium would be rewarded by putting the nice 
comparison of the DMRG results with a field-theoretical 
approach (which is a fundamental issue of Ref. I2T 
more solid grounds. 



2. On-resonance e = ±V/2 

Quite intuitively, the current through the resonant 
level vanishes if the latter is energetically placed above 
(or below) of cither one of the chemical potentials (e ^> 
V) while featuring the power-law decay described in the 
previous Section for small e <C V (see Fig. [5]). If the 
impurity position is fixed to the so-called resonance con- 
dition e = ±V/2 which determines the crossover between 
the aforementioned limits, the analytic expression (I5C 
takes the form 



J 



1 



1 



2U 



+ 0(U 2 ), (52) 



for voltages V Tk where the arc-tangent can be re- 
placed by its asymptotic value. However, even if V is 



orders of magnitude larger than the equilibrium energy 
scale Tx , the current does not necessarily exhibit a power 
law since the exponent a rcs becomes small for U 0. 
As before, we can verify this result numerically and ex- 
tract the quantity a ros beyond linear order by fitting 
to the form of Eq. (f52jl .— Within our (reservoir-cutoff) 
FRG approximation scheme, a Tes is equal to the off- 
resonance exponent aj (see the inset to Fig. E]) - Whereas 
those observations are altogether in complete agreement 
with recently-published real-time renormalization group 
results^ they clearly contradict the intuition that the 
voltage can always be interpreted in terms of an infrared 
cutoff associated with universal power laws^ 

While for vanishing impurity energies both functional 
RG cutoff schemes describe the same non-equilibrium 
physics, the violation of causality within the sharp cutoff 
approach in presence of finite e and voltages V leads to 
severe artifacts (e.g., deviations from power-law behav- 
ior for e <C V) already for fairly small Coulomb interac- 
tions. These violations originate from the diagonal self- 
energy components which do not flow at e = 0. Thus, 
the reservoir cutoff scheme is clearly superior in inves- 
tigating particle-hole asymmetric impurity positions for 
the IRLM beyond linear response. 



D. The left-right asymmetric IRLM 

So far, we have modeled the resonant level to be cou- 
pled symmetrically to both bath, and this scenario is 
generically employed in previous worksi 2 i 26 ' 34 i 35 How- 
ever, the FRG flow equations introduced in Sec. IIIII can 
be generalized straightforward for different Coulomb in- 
teractions Ul 7^ Ur and hopping amplitudes t' L 7^ f' R 
to the left and right side, respectively Following the 
same arguing as in the symmetric case, one can analyti- 
cally show that to leading order the current in the limit 
r > V > t' L Rl e takes the form 



J 



1 



'Tk 
, V 



'T K 



(53) 



where the linear-response low energy scale Tk and the 
asymmetry parameter c are given by Tk ~ ^'jJ-'r an d 
c = t' L /t' R to zcroth order. For Ul ^ Ur, J is thus not 
governed by a power law even if V 3> Tk is large, and the 
voltage can again not be interpreted as an infrared cutoff. 
A more detailed discussion of the two-channel interacting 
resonant level model with asymmetric couplings can be 
found in Ref. k32l 



VI. CONCLUSIONS 

In this work we have studied zero-temperature steady- 
state transport properties of the two-channel interact- 
ing resonant level model in presence of an arbitrary 
bias voltage V. Beyond linear response, the functional 
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rcnormalization group in Keldysh frequency space can be 
used to compute the self-energy associated with the lo- 
cal Coulomb interaction between the isolated level and 
the two baths of delocalized states. We truncate the 
infinite hierarchy of functional flow equations to lead- 
ing order, rendering the FRG an approximate method to 
calculate effective system parameters. Despite the simple 
(Hartree-Fock-like) nature of the resulting scheme, trans- 
port properties of the IRLM can be obtained to satisfying 
agreement with dynamical matrix renormalization group 
data both in and out of equilibrium. This allows for 
a thorough investigation of the special case where the 
bandwidth T of the leads becomes large. For voltages 
much smaller than the linear-response energy scale Tk, 
this so-called scaling limit is characterized by universal 
power-laws. E.g., the charge susceptibility is governed 
by x" 1 ~ £ /2 ~ Q x Tk, with t' being the local hop- 
ping to the leads. The exponent a x = 2U/irT + 0{U 2 ) 
can be computed numerically and analytically and agrees 
with prior results to leading order. Far from equilibrium 
(V 3> Tk), the current decays as a power-law J ~ V~ aj 
with the voltage if the impurity energy is small (e <C V) 
while featuring more complex behavior if the latter is 
pinned close to either one of the chemical potentials of the 
bath (e = ±V/2). Whereas in the former case the volt- 
age can be interpreted in terms of a simple infrared cut- 
off (which automatically appears within the FRG frame- 
work), the same does not hold close to the resonance 
condition e = ±V/2. 

From the methodical point of view, we have estab- 
lished the functional renormalization group as a sim- 
ple tool to compute effective (Hartree-Fock-like) param- 
eters incorporating aspects of non-equilibrium physics of 
quantum impurity systems. The latter particularly holds 
for a recently-proposed way of implementing an infrared 
cutoff in Keldysh frequency space, which can be inter- 
preted in terms of an additional reservoir whose cou- 
pling strength flows from infinity to zero and which does 
not suffer from symmetry violations specifically associ- 
ated with non-equilibrium (such as causality). In gen- 
eral, however, the simple approximation obtained from 
truncating the infinite hierarchy of FRG flow equations 
is limited to treat small to intermediate values of the 
Coulomb interaction only. Extending the method to the 
strong-coupling limit is subject to future work. 
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t' 2 /r 2 

FIG. 7: (Color online) Upper panel: The same as in Fig.0 but 
additionally showing FRG results obtained from the second 
order truncation scheme outlined in the Appendix. Lower 
panel: The logarithmic derivative dln(F /xt' 2 )/dln(t' 2 /Y 2 ) 
obtained from the FRG flow equation (J25J) as well as from a 
(frequency-independent) second-order generalization (see the 
main text for details) in comparison with numerical renor- 
malization group data. Both for the NRG as well as for the 
first order FRG scheme, this quantity is constant over orders 
of magnitude as the charge susceptibility \ i s governed by a 
power law of the bare hopping amplitude t' . 

Appendix: Second order functional RG 

In this Section, we briefly discuss how the functional 
RG scheme employed to derive the results of Sees. IIVI 
and [V] can be extended by accounting for the flow of the 
two-particle vertex. Whereas we focus solely on the case 
of linear response, a detailed presentation of second-order 
FRG data in the context of the Anderson impurity model 
out of equilibrium can be found in Ref. l2ll . 

Schematically, the flow equation of the two-particle 
vertex is given b y 18 ' 53 

9 A72 A ~ Tr S A 7 2 A G A 72 A + T ( 7 A ) , (54) 

and after neglecting the contribution of the three-particle 
vertex 7^ (i.e., truncating the infinite hierarchy to sec- 
ond order) one obtains a finite set of differential equa- 
tions which can in principle be solved numerically by in- 
troducing a discretization of Matsubara frequency space. 
For the single impurity Anderson model in equilibrium, 
this procedure (which is explained extensively in Ref. I20I ) 
leads to systematic improvements at small to interme- 
diate Coulomb interactions and additionally allows for 
computing energy-dependent properties such as the local 
density of states or the finite-temperature conductance^ 
but cannot describe strong-coupling physics such as the 
appearance of the exponentially small Kondo energy 
scaled 

For the problem at hand, one can pursue the same 
course of action and numerically solve the full (frequency- 
discretized) two-particle and self-energy flow equations 
(f54"f and (fTTj) . respectively. For the sake of numerical re- 
sources (and given the outcome) it is, however, desirable 
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to devise some simplifications. Here, we only account 
for density-density (nearest and next-nearest neighbor) 
interaction terms but exemplary ensured that addition- 
ally considering the flow of 'correlated hoppings' of the 
type d\d 2 d\d 1 docs not quantitatively alter our results 
for all cases of interest (particularly for the parameters 
of Fig. [3] as well as for the equilibrium exponent a x 
at [/ < r). Moreover, we approximate the frequency- 
dependence of the two-particle vertex by introducing 
three bosonic frequencies intrinsically connected to the 
three types of (particle-particle, particle-hole and hole- 
particle) diagrams appearing on the right-hand side of 
Eq. ([M]) and discard all but these specific frequency 
dependencies (this is again outlined in more detail in 
Ref.[20|). 

As for the Anderson model, taking into account the 
flow of 7^ in the so-specified way leads to systematic 
improvements of physical quantities such as the linear- 
response conductance (see the upper panel of Fig. [7]). 
In addition, Ward identities (connected, e.g., to differ- 
ent equivalent ways of computing the average occupa- 
tion number; see Ref. [20|) are in general violated by the 



FRG (which is a non-conserving approximation) but ful- 
filled to much greater accuracy for intermediate Coulomb 
interactions within the second-order scheme. In con- 
trast, the fundamental scaling-limit power-law behavior 
which is described by the differential equations (|25t - 
([27| in complete agreement with NRG reference data 
is no longer captured if the flow of 7^ is accounted for 
(see the lower panel of Fig. [7]). This already manifests if 
the frequency-dependence of 7^ is completely discarded 
(and the second-order flow equation describes effective 
nearest- and next-nearest neighbor Coulomb interactions 
only). This is again analogous to the Anderson model 
where an exponential energy scale contained within the 
most simple (Hartree-Fock-like) approximation scheme^ 
is no longer captured by the (more elaborate) second- 
order approach ! 20 ' 21 Thus, generalizing the FRG flow 
equations cannot be achieved by straightforward trun- 
cating at second instead of first order, and another strat- 
egy on how to tackle strong-coupling physics of quan- 
tum impurity systems using the functional renormaliza- 
tion group needs yet to be devised. 
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